Magnetic field induced quantum phase transition of the S = \/2 antiferromagnet 
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The magnetic properties of alkali- metal peroxychromate K2NaCr08 are governed by the S — 1/2 
pentavalent chromium cation, Cr^+. Specific heat, magnetocalorimetry, ac magnetic susceptibility, 
torque magnetometry, and inelastic neutron scattering data have been acquired over a wide range of 
temperature, down to 60 mK, and magnetic field, up to 18 T. The magnetic interactions are quasi- 
two-dimensional prior to long-range ordering, where Tat = 1.66 K in iif = 0. In the T ^ limit, 
the magnetic field tuned antiferromagnetic-ferromagnetic phase transition suggests a critical field 
He — 7.270 T and a critical exponent a — 0.481 ± 0.004. The neutron data indicate the magnetic 
interactions may extend over intra-planar nearest-neighbors and inter-planar next-nearest-neighbor 
spins. 

PACS numbers: 75.40.Cx, 75.30.Kz, 75.50.Ee 



The study of the similarities and differences between 
classical and quantum phase transitions has a rich his- 
tory, and recent research has focused on systems in re- 
duced dimensions ii^J^J^ For magnetic systems, the model 
Hamiltonian is often written as 
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where the first sum over is often restricted to 

nearest-neighbors (nn) but may sometimes extend to in- 
clude next-nearest neighbors (nnn), J > indicates an- 
tiferromagnetic interactions, and the Zeeman (second) 
term includes the presence of the externally applied mag- 
netic field H and an anisotropic Lande ^-tensor. It is 
important to emphasize that both the spin and spatial 
dimensions are crucial elements of any model, and in this 
context, the critical exponents can be considered finger- 
prints for classifying the nature of the transition. In par- 
ticular, the order parameter describing the transition can 
be characterized by (Tc — T)^ as the temperature T ap- 
proaches the critical temperature Tc or, in the T ^ 
limit, by {Hc — H)^ for a quantum phase transition tuned 
by the magnetic ^eldX^^d^^^^^^^^ 

The purpose of the present work was to consider the 
curious case of the pentavalent chromium cation, Cr^+, 
an 5 = 1/2 (3(i^) ion that is not a common choice for con- 
densed matter research. ^^^^^^^^^^^^^^^^'^^^^^^^^^^ Specifically, 
the alkali- metal peroxychromate K2NaCr08 crystallizes 
with orthorhombic symmetry (space group Phcm) with 
the lattice parameters a = 8.5883 A, 6 = 7.9825 A, and 
c = 9.2432 A at 173 K.^"^ The Cr^+ ion is bonded to 
four peroxy (02)^~ ligands in a dodecahedral geometry, 
and the ions form the three dimensional network shown 
in Fig. l^i^ii^ The data from EPR studies showed the 
^-tensor is anisotropic, with = 1.9851, gy = 1.9696, 




FIG. 1: (Color online) Spin network in K2NaCr08 as seen 
along the 6-axis, where only the Cr^^ ions are shown for clar- 
ity. The vertical planes are the mirror planes containing Ji 
(soUd lines) and the J2 (dashed lines) paths connecting the 
nn ions. Exchange paths corresponding to nnn inter- 

actions are labeled as J3 (short-dashed lines) and J4 (dotted 
lines) . 



and gx = 1.9636, the electronic ground state is d^^, 
and the magnetic z-axis is along the [110] crystal di- 
rection. During studies of the static magnetic suscep- 
tibility, Xdc{H = 0.3 T,1.8 K < T < 20 K), and the 
specific heat, C(0 < < 9 T, 1.8 K < T < 10 K), 
Cage and Dalal^^ observed ordering signatures in poly- 
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crystalline samples. More specifically, while a rounded 
maximum was observed in Xdc{H = 0.3 T, T ^ 2.29 K), 
the C{H 0,T ^ 1.8 K) data appeared to possess a 
A-like peak. These signatures indicate the presence of 
low dimensional antiferromagnetic behavior and the on- 
set of long-range antiferromagnetic ordering. More re- 
cently, Choi et al}^ performed low temperature magne- 
tization and ^^Na NMR measurements and identified a 
quantum phase transition from an antiferromagnetically 
ordered state to a ferromagnetically polarized state at 
critical field of about 7.4 T, with (3 = 0.44 ± 0.05 
in = 5 T and = 0.53 ± 0.05 in H = 6.8 T. Since 
mean- field theory yields /3 = 1/2 for magnetization, these 
workers concluded that the high field magnetization of 
K2NaCr08 is mean-field-like. 

Herein, the results of specific heat, magnetocaloric ef- 
fect, ac magnetic susceptibility, and torque magnetome- 
try are reported to construct the magnetic phase diagram 
of K2NaCr08, with particular emphasis on the T ^ 
limit of the magnetic field induced quantum phase tran- 
sition. The data indicate the magnetic interactions are 
quasi-two-dimensional prior to long-range ordering. In- 
elastic neutron scattering studies in zero magnetic field 
confirm the antiferromagnetic transition occurs without 
any low temperature structural transitions and suggest 
the strongest magnetic interactions may include both 
nn and nnn ions. Furthermore, in the critical regime, 
the specific heat and magnetocaloric effect systematically 
identify the phase transition, while the magnetic mea- 
surements, although consistent with each other, yield a 
slightly shifted trend as T ^ 0. 

Single crystals of K2NaCr08 were prepared by modi- 
fication of the method of Riesenfeld.^^ The synthesis in- 
volves the reduction of Cr^+ by 30% H2O2 at 5°C in the 
1:1 molar mixture of aqueous NaOH and KOH. Specific 
heat and magnetocaloric data were collected on single 
crystals with H applied along the c-axis. The specific 
heat studies employed the relaxation technique, 150 mK 
< T < 20 K, with < 18 T using two different instru- 
ments. The magnetic specific heat, Cmag, was obtained 
by subtracting the zero-field Debye contribution, namely 
Cphonon = (1.85 X 10-3) J/(mol K^), from the total 
specific heat at all fields, and was integrated to yield the 
magnetic entropy, SmagiT), Fig. 2. The magnetocaloric 
data, 150 mK < T < 600 mK, were collected in i7 < 18 T 
with a field sweep rate of 0.2 T/min, Fig. 3a. For the ac 
susceptibility, Xac, a sample consisting of a random ar- 
rangement of microcrystallites (8.44 mg) was mounted, 
using GE varnish dissolved in ethanol, on a silver strip 
and wrapped in Kapton tape. The measurements were 
performed in a homemade dilution refrigerator probe us- 
ing a standard set of inductively coupled coils operating 
at 1.114 kHz, Fig. 3b. Torque magnetometry, performed 
down to 280 mK and up to 18 T, utilized a single crys- 
tal sample (1.2 mg) mounted using Stycast epoxy on a 
0.005 inch thick CuBe cantilever soldered to a post. The 
torque experienced by the sample deflects the cantilever, 
and the deflection was measured capacitively through a 
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FIG. 2: (Color online) The temperature dependence of 
the magnetic entropy, 5'mag(T), at several representative 
magnetic fields. The inset shows the temperature depen- 
dence of the magnetic specific heat divided by temperature, 
Cmag(T)/T, at the same magnetic fields. Cmag(T) was ob- 
tained by subtracting a lattice contribution from the total 
specific heat, see text. Smag{T) was obtained by integrating 
the data shown in the inset. As T increases, Smag{T) slowly 
approaches the R\n2 value, designated by the dashed line, as 
expected for a 5* = 1/2 system. In iif = 0, about 60% of the 
entropy is acquired above Tn , indicating significant contribu- 
tions from short range correlations. 



manual capacitance bridge excited at 30 V and 5 kHz, 
Fig. 3c. Zero field neutron scattering work down to 1.4 K 
used a powder-like sample, generated by grinding micro- 
crystals, that was mounted on the NIST-CNS disk chop- 
per spectrometer (DCS) tuned to a wavelength of 4.8 A. 

The distinct A anomaly observed in the magnetic spe- 
cific heat (Fig. 2) at Tat = 1.66 K and = in- 
dicates that K2NaCr08 undergoes long-range antiferro- 
magnetic ordering (LRAFO) as conjectured from earlier 
specific heat measurements.^^ The LRAFO is suppressed 
in the presence of a magnetic field, as Tn shifts slowly 
to lower temperatures until H reaches about 5 T and 
rapidly thereafter until the LRAFO is completely sup- 
pressed around 7.3 T. For > 5 T, the Schottky-type 
anomaly overlapping the LRAFO at zero field becomes 
the dominant feature of Cmag(^), and the corresponding 
position in temperature of the heat capacity maximum 
shifts to higher temperatures. Consequently, the posi- 
tion of the phase transition is harder to identify. On the 
other hand, the magnetocaloric effect produces peaks in 
the up and down sweep directions, and these features 
correspond to the magnetic phase transition, while the 
center of these anomalies is identified as the transition 
point. Fig. 3. In Fig. 4, the transition points identified 
by these two methods are shown to yield a single trend 
of the phase boundary between the LRAFO and the po- 
larized ferromagnetic state. 
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FIG. 3: (Color online) (a) The magnetocaloric effect curves 
for up and down field sweeps near 230 mK. (b) The ac sus- 
ceptibility, Xac{H) at 69 mK and 370 mK. (c) The isothermal 
magnetic torque measured as a function of H. In each panel, 
the arrows designate the positions taken as the field induced 
phase transition. 



With respect to the magnetic studies using Xac and 
torque, the signatures designating the critical field of the 
transition are sharpest as T ^ 0. For Xac, the inflec- 
tion point between the low and high field regions is often 
used, Fig. 3b. Alternatively, these data are integrated to 
yield the isothermal field dependence of the magnetiza- 
tion, M(iJ), but an extrapolation protocol must be used 
as the location of the critical field is often smeared*^*^*^ 
In our case, both routines provide essentially the same 
results that are plotted in Fig. 4. For the torque mea- 
surement, the peak of the response appears on a smooth, 
nearly quadratic background. Fig. 3c, and is taken as the 
location of the transition. Fig. 4. It is interesting to note 
that both techniques indicate a phase boundary that is 
slightly shifted from the one given by the specific heat 
and magnetocaloric effect results. Fig. 4. 

Since the response of the specific heat unambiguously 
marks the transition and the magnetocaloric effect data 
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FIG. 4: (Color online) Magnetic phase diagram of K2NaCr08. 
The solid line corresponds to Eq. 2 with He = 7.27 T and 
a = 0.481 ± 0.004. The right panel provides an expanded 
view near Hc{T 0), where the specific heat and magne- 
tocaloric transition points are systematically different than 
the ones extracted from standard analysis of the magnetic 
measurements, Xac and torque. 



identify the same critical field, these results were used to 
fit to 



(2) 



Care must be used when fitting Eq. 2, as both He 
and a must be resolved i^i^ii^iii Consequently, the "T- 
window" analysis and empirical convergence methods^' 
were employed, by systematically varying the maximum 
temperature of the fits from 1 K to 230 mK, to extract 
the values of He = 7.27 T and a = 0.481 ± 0.004. In 
detail, to extract i^c, the low temperature (T < 1 K) 
phase boundary was fit to Eq. 2 for various a values in a 
temperature window, Tw^ of decreasing size. For exam- 
ple, the phase boundary with a window size of 20 data 
points was fit to Eq. 2 with He and A as fitting param- 
eters, while holding a = 0.48. Then the window was 
successively reduced by removing the highest tempera- 
ture point, and the fitting procedure was repeated. The 
resulting He values can then be plotted as a function of 
the upper limit of Tw^ Fig. 5. The whole fitting pro- 
cedure is then repeated for various fixed a values. The 
value of He = 7.2696 ± 0.0002 T was then obtained by 
taking a weighted average of the linearly extrapolated 
HeiTw 0) values corresponding to various a values. 
This convergence of He to a single value allows a more 
precise determination of a than possible with a three pa- 
rameter fit in which He, a, and A are simultaneous fitting 
parameters. The critical exponent a is then determined 
by following the same procedure, where the T-window 
analysis is performed when fixing He = 7.2696 T, Fig. 5, 
and the result is a = 0.481 ±0.004. Establishing the crit- 
ical exponent of the transition begs the questions about 
the details of the network of magnetic interactions, so 
inelastic neutron scattering experiments were performed. 
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FIG. 5: (Color online) Estimates of (top panel) the critical 
field He and (bottom panel) the critical exponent a obtained 
from the "T-window" analysis of the transition, see text for 
detailed description. The results are He = 7.2696 =b 0.0002 T 
and a = 0.481 ± 0.004. 
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FIG. 6: (Color online) Inelastic neutron scattering spectra at 
(a) 1.725 K and (b) 1.7 K. 



The neutron scattering pattern at 1.725 K (Fig. 6a) 
reveals a magnetic excitation above = 0.2 meV. 
After decreasing the temperature to 1.7 K (Fig. 6b), 
this excitation becomes sharper with a well-defined mode 
near = 0.5 meV. When the magnetic excitations are 
summed over all Q, the broad peak near /S.E = 0.5 meV 
sharpens for T < 1.725 K. The temperature dependence 
of the integrated area of the peak is consistent with the 
LRAFO observed in the zero-field specific heat. 

Another noteworthy feature is that the excitation be- 
low the transition shows clear Q dependence, which has 
higher intensity at low Q. This magnetic contribution 
is consistent with an expression used to model the scat- 
tering from powders, where short-range, isotropic spin 
correlations persist to the first few coordination shells of 
neighboring spins, namely 

HQ) - l^{s,s,)^-- ^ , (3) 

where the final result is for the case of spins correlated 
over the nearest neighbors. Consequently, only one value 
R of Vij (the distance between spins at sites i and j) is em- 
ployed in our fitting procedure, where Eq. 3, multiplied 
by the magnetic form factor of the Cr^+ ion, has been 



used to fit to the net intensity shown in Fig. 7a. This 
analysis yields = 6 A, which is close to the nn and 
nnn Cr- • • Cr distances that are in the range 5.8-6.2 A. 
Ergo, the magnetic correlations are short range, extend- 
ing over nn and nnn neighbors. Fig. 1. Furthermore, 
with decreasing temperature, the lattice Bragg peaks do 
not show obvious splitting or intensity change, thereby 
excluding the possibility of a structural change. On the 
other hand, new magnetic Bragg peaks with low inten- 
sity develop below 1.725 K (Fig. 7b) and are consistent 
with the ordering of a small moment at ^ 1.7 K. 

Equipped with the result that short range magnetic in- 
teractions appear to extend over nn and nnn ions, the un- 
derlying network of these interactions can be reexamined. 
Previous theoretical and experimental studies established 
that although most of the unpaired electron density is 
on the Cr^+ ion, a significant portion of it can be found 
on the bound (62)^" ions.— i^i^*^ Consequently, the 
potential superexchange pathways between Cr^+ ions in 
K2NaCr08 can be expected to involve the (62)^", Na+ 
and ions. Since the bond angle is an important factor 
in determining the ferromagnetic versus antiferromag- 
netic nature of the interactions,— ^^^^ the paths with 
larger O- • • Na/K- • • O angle are selected when more than 
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FIG. 7: (Color online) (a) A cut of the magnetic excitation 
near = 0.5 meV at T = 1.4 K. The Q dependence of the 
intensity (open squares) is fit to Eq. 3 (solid line), (b) Ap- 
pearance of the magnetic peaks at 1.7 K, after the data at 
2 K have been subtracted as the lattice standard. 

one possible exchange pathway for a given Cr^+- • • Cr^+ 
distance are present. This analysis leads to the results 
given in Table I, which tabulates the most influential ex- 
change interactions between nn (the intra-planar Ji and 
J2)i^ and nnn (the inter-planar J3 and J4) Cr^+ ions, 
Fig. 1. While the Ji path involves Na+ ions with a 
Na- • • O distance of ~ 2.4 A and a bond angle ZONaO of 
142°, the J2 path involves K+ ions with K- • • O distance 
of - 3.0 A and an angle ZOKO of 166°. The J3 path 
involves Na+ ions with a Na- • • O distance of ~ 2.4 A and 
a bond angle ZONaO of 166° and the J4 path involves 
K+ ions with a K- • • O distance of ~ 2.6 A and a bond 
angle ZOKO of 131°. All of the data suggest the intra- 
layer interactions are somewhat stronger than the inter- 
layer ones, since the two-dimensional correlations de- 
velop before three-dimensional long-range order appears. 
Assuming a square lattice model, the intra-planar ex- 
change interaction, Jintra = ('^i + '>^2)/2, can be estimated 



from the data. For example, fitting Cmag{T^H = 18 T) 
yields Jintra ~ 2.1 K, whereas He = 4Jintra/^MB, gives 
Jintra ^ 2.4 K. Both of thcsc estimates are close to the 
value of 2.44 K obtained by fitting Xdc{T) using a two- 
dimensional high temperature series expansion. 

In summary, the low temperature, high magnetic field 
phase diagram of K2NaCr08 tias been established. In the 
T ^ limit, the magnetic field tuned antiferromagnetic- 
ferromagnetic phase transition suggests a critical field 
He = 7.270 T and a critical exponent a = 0.481 ± 0.004. 
The value of this critical exponent, Eq. 2, depends on 
the universality class of the quantum phase transition, 
where a = 2/3 is expected for Bose-Einstein conden- 
sation (BEC)^^^'^^^ and a = 1/2 is predicted for Ising- 
like spins The present data do not allow us to 
unambiguously differentiate between these two possible 
descriptions, since a broad range of exponents is ex- 
pected if the transition is not probed at sufficiently low 
temperatures.^ The neutron data suggest the magnetic 
interactions extend over intra-planar nearest-neighbors 
and inter-planar next-nearest-neighbor spins. For a 
deeper understanding of these results, especially the 
value of the critical exponent, inelastic neutron scattering 
experiments should be performed on single crystals that 
are cooled to millikelvin temperatures in the presence of 
the magnetic field capable of spanning He. Fluctuation- 
induced heat release^ and NMR^i studies are alterna- 
tives that might provide decisive information. 

All experimental studies, except the ac susceptibil- 
ity and neutron measurements, were performed at the 
National High Magnetic Field Laboratory (NHMFL), 
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Agreement Grant No. DMR-0654118 and by the State 
of Florida. We thank C. D. Batista, K. Ingersent, and 
D. L. Maslov for insightful conversations or communica- 
tions, and we are grateful to Y. Qiu and J. R. D. Cop- 
ley for help with the neutron measurements. This work 
was supported by NSF, in part, via DMR-0520481 (SN), 
DMR-0701400 (MWM), and DMR-0803516 (YL). 



TABLE I: Possible magnetic exchange pathways between 
Cr^^ in K2NaCr08. The labels for the atoms is the same 
as used by Cage and Dalal^^ and the Ji and J2 designations 
are the same as used by Choi et al}^ 

Label Cr- • • Cr (A) Pathway 

5^837 01-Na2-02 

J2 5.888 02-K1-02 

J3 5.974 06- • • Na2- • • 06 

J4 6.243 03-K1-.-03 
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